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Abstract 



^ \ Hyperfine quenching rates of the lowest-energy metastable ^Pq and states of Mg, Ca, Sr, and 

W ■ Yb atoms are computed. The calculations are carried out using ab initio relativistic many-body 

■ methods. The computed lifetimes may be useful for designing novel ultra-precise optical clocks 

and trapping experiments with the fermionic isotopes. The resulting natural widths of the 
^ ^! -^'S'o clock transition are 0.44 mHz for ^^Mg, 2.2 mHz for ^^Ca, 7.6 mHz for '^'^Sr, 43.5 mHz for 

^''^Yb, and 38.5 mHz for ^'^^Yb. Compared to the bosonic isotopes, the lifetime of the states in 
fermionic isotopes is noticeably shortened by the hyperfine quenching but still remains long enough 
for trapping experiments. 

^ ■ PACS numbers: 32.70.Cs, 32.10.Fn, 31.15.Ar 
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I. INTRODUCTION 



This work is motivated by emerging experiments with cold divalent atoms Mg, Ca, Sr 
and Yb 0. For example, the recently attained Bose-Einstein condensate of the ground-state 
Yb may offer new insights into the physics of degenerate quantum gases due to a vast 
number of available isotopes and relative simplicity of molecular potentials. As to the 
metastable states (see Fig.^, it was realized that the non-scalar nature of the states may 
be used to overcome the unfeasibility of magnetic trapping of the spherically-symmetric ^5*0 
ground states Knowing radiative lifetimes of the other, ^Pq, metastable states is 

required in developing the next generation of ultraprecise optical atomic clocks 0, H S ■ 
Here the clockwork is based on cold atoms confined to sites of an engineered optical lattice. 
The lifetime determines the natural width of the clock transition between the ground and 
the ^Pq state. 

'P. (nsnp) ■ 



(nsnp) 



'Sn(ns2) 



FIG. 1: Lowest-lying energy levels of Mg (n=3), Ca (n=4), Sr (n=5), and Yb (n=6), relevant to 
the radiative decay of the nsnp^Po,2 states. The hyperfine quenching predominantly is caused by 
the admixture of the nsnp^Pi and nsnp^Pi states. 

For all bosonic isotopes of Mg, Ca, Sr, and Yb, the nuclear spin I vanishes and these 
isotopes lack hyperfine structure. For bosonic isotopes the ^Pq state may decay only via very 
weak multi-photon (e.g., El-Ml) transitions. However, for fermionic isotopes (Table H]), 
/ 7^ 0, a new radiative decay channel becomes available due to the hyperfine interaction 
(HFI). The HFI, although small, admixes atomic levels of the total angular momentum 
J = 1 thus opening an electric-dipole branch to the ground state. The resulting HFI- 
induced El decays do determine the lifetimes of the ^Pq states. As to the ^P2 states, here 
the single-photon decays are allowed, but being of non-El character, are very weak. The 
lifetimes are long and range from 15 seconds for Yb to 2 hours for Ca 11, 2|. As we 



demonstrate here, depending on an isotope, the hyperfine quenching of the P2 states is 
either comparable to or is much faster than the small non-El rates. 

A detailed theoretical analysis of the hyperfine quenching has been limited so far to 
astrophysically-important He, Be and Mg and their isoelectronic sequences |13„. .14. 15, 16 J . 
The hyperfine quenching of the ^Pq states of Sr and Yb has been estimated in Refs. [7|, |lO| . 
Here we carry out ab initio relativistic many-body atomic structure calculations to extend 
and refine these previous studies. We find that the resulting natural widths of the ^Pq —^Sq 
clock transition are 0.44 mHz for ^^Mg, 2.2 mHz for ^'^Ca, 7.6 mHz for ^^Sr, 43.5 mHz 
for ^^^Yb, and 38.5 mHz for ^''^Yb. Compared to the bosonic isotopes, the lifetime of the 
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TABLE I: Nuclear parameters of the stable fermionic isotopes of Mg, Ca, Sr, and Yb. Here / are 
the nuclear spins and n/fJ-N are the nuclear magnetic moments expressed in units of the nuclear 
magneton njy. 



Isotope / fJ-ffJ-N 

25 Mg 5/2 -0.85546 

^^Ca 7/2 -1.31727 

^^Sr 9/2 -1.09283 

I'^^Yb 1/2 0.4919 

I'^^Yb 5/2 -0.6776 



states in fermionic isotopes is noticeably shortened by the hyperfine quenching but still 
remains long enough for trapping experiments. 

The paper is organized as follows. First, in Section |n] we derive the hyperfine quenching 
rates using perturbation theory. The solution of many-body atomic problem and numerical 
details are given in Section ITTTl Finally, we present the results, compare with the previous 
calculations, and draw the conclusions in Section HV1 Unless noted otherwise, atomic units 
[h = \e\ = rrie = 1) are used throughout. 



II. DERIVATION OF HYPERFINE QUENCHING RATES 

In the presence of nuclear moments, the total electronic angular momentum J no longer 
remains a good quantum number. The atomic energy levels are characterized instead by the 
total angular momentum F = J + I. Nevertheless, the coupling of the electronic and the 
nuclear momenta is small and in this section we employ the first-order perturbation theory 
in the magnetic-dipole hyperfine interaction to compute the modified atomic wave functions 
of the levels. With these perturbed wave functions the hyperfine quenching rates are 
obtained with the conventional Fermi golden rule. 

Before proceeding with the outlined derivation, we notice that in this problem there are 
two types of perturbations: the hyperfine interaction and the interaction with the electro- 
magnetic field. Here we treat the HFI as the dominant interaction, determine the hyperfine 
structure first, and as the next step compute the lifetimes. This approach is valid as long as 
the radiative width of the ^Pi level is much smaller than the fine-structure intervals between 
the components of the ^Pj multiplet 3|. We verified that this inequality holds for all the 



atoms under consideration. 

We develop the formalism in terms of the hyperfine states \'~f{IJ)FMp). Here the angular 
momenta I and J are conventionally coupled to produce a state of definite total momentum 
F and its projection Mp, and 7 encapsulates all other atomic quantum numbers. In the 
first order of perturbation theory in the hyperfine interaction, Hupi, the correction to the 
hyperfine sub-level \-y{IJ)FMF) of the metastable state I7J) reads 



7' J' 



where E (7 J) are the energies of atomic states. In the above expression, we have taken into 
account that -ffHFi is a scalar, so the total angular momentum F and its projection Mp are 
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conserved. In general, the hyperfine coupling Hamiltonian, Hupi, may be represented as a 
sum over multipole nuclear moments Ai^'^^ of rank k combined with the even-parity electronic 
coupling operators T^^^ of the same rank so that the total interaction is rotationally and P- 
invariant. For the hyperfine quenching of the J = states via couplings to rapidly decaying 
J = 1 states, we may truncate the Hupi at the magnetic-dipole part 

H,,,= (t!l.T^A. (2) 

Here /i^^-* is the operator of the nuclear magnetic moment, fijy is the nuclear magneton, and 
is a relevant operator acting in the electronic space. 

While for the nsnp^Po the truncation of the HFI Hamiltonian at the nuclear magnetic- 
dipole contribution ^ is rigorously justified due to the selection rules, for the nsnp^P2 
states such a truncation requires a special consideration. Indeed, the next leading term 
in the multipole expansion of H^pi is due to the nuclear electric quadrupole moment. The 
associated electronic tensor of rank 2 can also admix levels of J = 1 symmetry and contribute 
to the hyperfine quenching. This quadrupole contribution does not vanish for isotopes with 
/ > 1, i.e., it may modify the nsnp^P2 for all the isotopes listed in the Table HJ except for 
^^^Yb. Our truncation of the HFI Hamiltonian will lead to a small "systematic" error for 
the nsnp hyperfine quenching rates and we will return to this point in the conclusion. 

With the magnetic-dipole contribution, Eq. Q, the required mixing matrix element in 
Eq. (HD is 

{^'{IJ')F'M'p\HufMIJ)FMf) = 6ff'Sm^m'^ 



^ (-1)-'- AyM±M±i) ^, ^} (y,'||r(..|M, (3) 

where n are the nuclear magnetic moments, compiled in Table EJ 

Given the correction to the wavefunction, Eq.(0), we derive the hyperfine quenching rate 
using the standard formalism. The rate of spontaneous emission (a — > b) for an electric- 
dipole radiation is 

Aa^, = ^u;l,\{a\T>\b)\\ (4) 

where a ~ 1/137 is the fine structure constant, Uab = Ea — Eb is the transition frequency, 
and D is the electric-dipole operator. Summing over all possible Fb and magnetic quantum 
numbers Mb of the final state, while disregarding small F-dependent energy correction, one 
obtains 

T 2Fa + 



For the case at hand, the initial state is the HFI-perturbed nsnp^Pj, [J ^ 1) state, and the 
final state is the ground ns^ ^S^ state. Taking into account Eq.(0), we arrive at the hyperfine 
quenching rate 



A ( 3p p 21c ^ 3//i V (2/+1)(J+1) f//l 
Ahfi \nsnp Pj; F ^ ns Sq) = -^^j ( — 1 j \j \ F] ' ' 
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with ujj = E (nsnp^Pj) — E {ns^ ^Sq) and the sums Sj defined as 

*Pj||T«||yj')(7'j'||D||; 

E J') - E {nsnp^Pj) 



{nsnp'^Pj\\T^^^\iJ'){iJ'\\D\\ns^^So) 



7',J' 



Notice that due to the electric-dipole selection rules, J' = 1. Also while the rate Ahfi 
depends on the nuclear parameters and the value of F, the sums Sj do not. In a particular 
case of the nsnp ^Pq states, the rate formula (jH)) may be simplified to 



Ahfi {nsnp 'Po; F = I ^ us' X) = ^u', (j^^ ^ l^o^ ■ (8) 

In the following section we describe the ab initio relativistic many-body calculations of the 
derived hyperfine quenching rates. 



III. SOLVING ATOMIC MANY-BODY PROBLEM 

The ab initio relativistic atomic-structure calculations employed here are similar to com- 
putations of electric-dipole amplitudes for the alkaline-earth ato ms |l7l| and hyperfine struc- 
ture constants and electric-dipole amplitudes for ytterbium [ij, Il8l|. Here we only briefiy 
recap the main features of this method. We consider Mg, Ca, Sr and Yb as atoms with two 
valence electrons outside the closed-shell cores. Strong repulsion between the two valence 
electrons is treated non-perturbatevely using the configuration-interaction (CI) method. The 
core-valence and core-core correlations are taken into account with the help of the many- 
body perturbation theory (MBPT) method. In the following we refer to this combined 
approach as the CI+MBPT method [l9| . 

In the CI-I-MBPT approach, the energies and the wave functions are determined from 
the eigenvalue equation in the model space of the valence electrons 

H,s{Ep)\%)=Ep\%), (9) 

where the effective Hamiltonian is defined as 

H,s{E) = Hpc + ^{E). (10) 

Here Hpc is the relativistic two-electron Hamiltonian in the frozen core approximation and 
E(i?) is the energy- dependent core-polarization correction. The all-order operator S(i?) 
completely accounts for the second order correlation correction to the energies. The omitted 
diagrams in higher orders may be accounted for indirectly by adjusting the effective Hamil- 
tonian 



17|,|20|. Namely, one introduces an energy shift S and replaces S(£') with Ti{E — 6). 
The parameter S is determined semi-empirically from a fit of the resulting theoretical energy 
levels to experimental spectrum. 

Using the effective Hamiltonian we find the wave functions of the ground and the ^Pj 
states. Further we apply the technique of effective all-order ("dressed") operators to cal- 
culations of the matrix elements. Technically, we employ the random-phase approximation 
(RPA). The RPA sequence of diagrams describes a shielding of externally applied field by 
the core electrons. This is the level of approximation employed here for electric-dipole ma- 
trix elements. The hyperfine, T^^\ matrix elements required more sophisticated approach: 
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TABLE II: Magnetic-dipole hyperfine structure constants A for the nsnp^Pi and nsnp^P2 states. 
The computed values are compared with the experimental data. 







A(3p°) (MHz) 


A{^P°) (MHz) 


25Mg 


This work 


-146.1 


-129.7 




Experiment 


-144.977(5) " 


-128.445(5) " 


43Ca 


This work 


-199.2 


-173.1 




Experiment 


-198.890(1) 


-171.962(2) " 


s^Sr 


This work 


-258.7 


-211.4 




Experiment 


-260.083(5) 


-212.765(1) 


171 Yb 


This work 


3964 


2704 




Experiment 


3957.97(47) 


2677.6 1 


173Yb 


This work 


-1092 


-745 




Experiment 


-1094.20(60) ^ 


-737.7 f 



Turio |2l|, ''Arnold et al. 0, ''Grundevik et al. 0, ''Heider and Brink"j2^], "Clark et al. 0, 

^Budick and Snir 0. 



for this operator we additionally incorporated smaller corrections (e.g., normalization and 
structural radiation; the details can be found in Ref. [Ts']). For the heaviest and more 
computationally demanding Yb, the corrections to the effective hyperfine operator tend to 
cancel y^], and we have simplified the calculations for Yb by using the bare T*^^) operator. 

To demonstrate the quality of the constructed wave functions and the accuracy of the 
effective-operator approach, in Table ITTl we present the calculated magnetic-dipole hyperfine 
structure constants A for the ^Pi_2 states. These constants are expressed in terms of ex- 
pectation values of Hn-pi- As seen from the Table ITTl the differences between the calculated 
and the experimental values, even for heavy Yb, do not exceed 1%. Further the sums Sj, 
Eq. ((7j), are computed in the framework of Sternheimer-Dalgarno-Lewis method 
the heart of this method is the recasting of the sums Sj in the form 

Sj = {nsnp^Pj\\T^^^\\5^), (11) 

where |(5\E') satisfies the inhomogeneous Schrodinger equation 

(i/eff - Einsnp'^Pj)) \6^>) = D,^\ns^ %) . (12) 

It is worth noting that because the effective operators act in the valence model space, the 
solution encompasses only the excitations of the valence electrons to higher valence 
states. The unaccounted for core excitations involve large energy denominators and we 
disregard their contributions. 



IV. RESULTS AND CONCLUSIONS 

To reiterate the discussion of the previous section, we carry out the calculations in several 
logical steps. First, we solve the CI+MBPT eigenvalue problem (jUJ and determine the 
ground and the nsnp^Pj state wavefunctions and energies. At the next step, we compute the 
dressed El-operator D^s and solve the inhomogeneous equation (fT^ . Finally, we calculate 
the required sums Sj, Eqs.ff7|l.(fTT|l. and determine the hyperfine quenching rates, Eq. (jH)). 
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The computed values of the isotope-independent sums 5*0 and 5*2 for Mg, Ca, Sr, and Yb 
are presented in Table IIIII The sums grow larger for heavier atoms due to increasing matrix 
elements of the hyperfine interaction (see Table ITT]) . A direct investigation of the sums shows 
that the contributions of both nsnp^Pi and nsnp^Pi intermediate states are comparable. 
The triplet state is separated by just a fine-structure interval from the metastable states, 
but its El matrix element with the singlet ground state vanishes non-relativistically. For the 
singlet state, the situation is reversed: compared to the triplet contribution, the involved 
energy denominator is much larger, but the electric-dipole matrix element is allowed. 



TABLE III: Sums Sj for the metastable ^Pq and ^P2 states. The values are given in atomic units. 





So,nsnp^Po 


S2 , nsnp ^P2 


Mg 


1.36 xlQ-^ 


2.30x10-5 


Ca 


3.56 xl0~5 


5.13x10-5 


Sr 


8.86 xl0~5 


1.27x10-^ 


Yb 


2.27 xlO""^ 


3.83x10-^ 



With the determined values of Sj and Eq.(jni), we obtain the hyperfine quenching rates for 
the metastable ^Pq and ^P2 states. The resulting rates are listed in Table ITVl The tabulated 
decay rates for the ^P2 states require some explanation. First of all, as follows from Eq.(jH)) 
the quenching rates depend on the total angular momentum F of the hyperfine substate. 
Although, in general, the total angular momentum F ranges from \ J — I\ to J + I, the 
6j-symbol in Eq.® imposes a stronger restriction, |/ — 1|<F<J-|-1. This requirement 
can be tracked to the selection rule for the electric-dipole transition amplitude between the 
ground state {Jg = 0,Fg = I) and the intermediate state which has the same F as the 
original hyperfine state (see Eq. ([[])). Keeping this restriction in mind, in Table ITVl we have 
listed the quenching rates only for such El-allowed values of F. 

We also remind the reader that in our analysis we have disregarded the contributions of 
the quadrupole and higher-order nuclear magnetic moments. While for the J = states 
this truncation is rigorously justified, for the J = 2 states the quadrupole contribution is 
generally present and becomes increasingly important for heavier atoms. Its relative role 
may be roughly estimated by forming a ratio of the relevant electric-quadrupole, B, and the 
magnetic-dipole. A, hyperfine-structure constants. The ratio B/A for the nsnp^P2 states 
is less or in the order of 0.1 for isotopes of Mg, Ca, and Sr, but is larger than unity for 
i73Yi3_ We expect that the quadrupole correction will be significant for ^''^Yb. However, for 
the ^^^Yb isotope, 7=1/2 and the nuclear moments beyond the magnetic-dipole moment 
vanish, justifying the validity of the truncation. 

Based on better than 1% accuracy of the ab initio hyperfine constants (Table |n} and 



energy levels [17|, |18[ we expect that the computed hyperfine quenching rates for the Pq 
states are accurate within at least a few per cent. For the ^P2 states of alkaline-earth 
atoms the main source of uncertainty is due to the neglected nuclear quadrupole moment 
contributions; the overall accuracy should be worse than that. 

In Table ITVl we also compare the computed rates with the results from the literature. For 
Mg the hyperfine quenching rates for the ^P2 state were estimated more than four decades 
ago by Garstang [13|. Our results are in a reasonable agreement with his values. Certainly, 
our calculations based on the modern ab initio relativistic many-bod y te chniques are more 
complete. For instance, in the calculation of the sums Sj, Garstang [l3[ kept only the two 
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TABLE IV: The hyperfine El-quenching rates for the metastable ^Pq and '^^2 states m sec~^. The 
rates depend on the total angular momentum F. The rates are compared with values by other 
authors, where available. 



Atom 


Transition rate 


F 


This work 


Other 


25]VIg 


AuFsi^Po - 


> 'So) 


5/2 


4.44 X 10"^ 


4.2 X IQ-^ " 




AnFs{^P2 - 




3/2 


2.31 X 10"^ 


1.4 X 10-"^ " 








5/2 


4.69 X 10"^ 


2.9 X 10-^ 








7/2 


4.95 X 10"^ 


3.1 X 10-^ 


43Ca 




> ^^o) 


7/2 


2.22 X 10^3 






^HFs(^i^2 - 




5/2 


1.02 X 10^3 










7/2 


1.81 X 10~3 










9/2 


1.74 X 10~^ 




s^Sr 






9/2 


7.58 X 10^3 


6.3 X 10-3 




AnFs{'^P2 - 


> ^^o) 


7/2 


4.13 X 10^3 










9/2 


6.86 X 10^3 










11/2 


6.27 X 10-3 




171 Yb 


AhfsCPo - 


> 'So) 


1/2 


4.35 X 10-2 


5.0 X 10-2 " 




AnFsi'^P2 - 




3/2 


9.18 X 10-2 




173Yb 


AKFsi^Po - 




5/2 


3.85 X 10-2 


4.3 X 10-2 




'^Garstang ^13.] . 


''Katori et al. J7], 


^Porsev et al. 





lowest-energy intermediate states ^Pi and ^Pi. This author has also employed the following 
El-matrix elements \(^So\\D\\^Pi)\ = 0.0058 a.u. and \{^So\\D\\^Pi)\ = 3.46 a.u., which 
are smaller than more accurate values [l^ of 0.0064(7) a.u. and 4.03(2) a.u., employed 
here. Our results for ^''Sr are in fair agreement with the estimate of Ref. Previously, 
we have estimated the quenching rates for Yb isotopes [l3 by summing only over the two 
lowest-energy excited states; the present result should be considered as more accurate. 

In Table IVl the calculated hyperfine quenching rates (maximum over hyperfine manifold) 
for the states are compared with the conventional electromagnetic transition rates. For 
Mg, Ca, and Sr these rates were calculated in Ref. [l2] and are due to Ml, M2, E2, and 
E3 multipole transitions. If the hyperfine quenching is allowed for a particular value of F, 
both rates contribute at a comparable level for Mg. For Ca and heavier atoms the hyperfine 
quenching becomes the dominant decay branch and determines the lifetime of the fermionic 
isotopes. 

It is worth mentioning that Yasuda and Katori [i^ have experimentally demonstrated 
that the blackbody radiation at 300 K quenches the ^Pj metastable states of Sr, significantly 
shortening their lifetimes. We verified that for other atoms, Mg, Ca, and Yb at T < 300 
K the blackbody radiation does not affect the lifetimes of the ^Pj states. We would like to 
emphasize, that the results tabulated in this paper are for the radiative decay rates due to 
the vacuum fluctuations of the electromagnetic field, i.e., for the ambient temperature of 
T = 0. The additional quenching by the temperature-dependent blackbody radiation [i^ 
should be also included in the total rate, especially for Sr isotopes. 

To summarize, here we employed the relativistic many-body methods to evaluate the 
hyperfine quenching rates for the metastable ^Pq and states of Mg, Ca, Sr, and Yb. The 
tabulated rates may be useful for designing novel ultra-precise optical clocks and trapping 
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TABLE V: Comparison of the hyperfine quenching rates (maximum over the hyperfine manifold 
in Table HVj) with the non-dipole rates for the states. The rates are given in sec~^. 



Atom 


hyperfine rate, max 


non-El rate " 


25]VIg 


4.95 xlO-"^ 


4.42 xlO-'' 


^3Ca 


18.1 xlO""^ 


1.41 xlO""^ 


«7Sr 


68.6 xlO"'' 


9.55 xlO"'' 


171 Yb 


9.18 xlO-2 


6.7 xlO-2 



"For Mg. Ca. and Sr the rates are from Ref. |12|| and for Yb from Ref. [111]. 



experiments with fermionic isotopes of metastable alkaline-earth atoms and Yb. The result- 
ing natural widths of the ^Pq —^Sq clock transition are 0.44 mHz for ^^Mg, 2.2 mHz for ^^Ca, 
7.6 mHz for ®^Sr, 43.5 mHz for ^^^Yb, and 38.5 mHz for ^^^Yb. Compared to the bosonic 
isotopes, the lifetime of the states in fermionic isotopes is noticeably shortened by the 
hyperfine quenching but still remains long enough for trapping experiments. 
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